Preprint typeset in JHEP style. - HYPER VERSION 



CERN-TH 2000-363 
IFIC/00-77 



Large lepton mixing and supernova 1987A 



M. KachelrieB 

TH Division, CERN, Geneva 

R. Tomas and J.W.F. Valle 

Instituto de Fisica Corpuscular - C.S.I.C. - Universitat de Valencia 

Ed. de Institutos de Paterna - Apartado de Correos 22085 - 46071 Valencia, 

Spain 



Abstract: We reconsider the impact of u e <->■ neutrino oscillations on the ob- 
served i> e signal of supernova SN 1987A. Performing a maximum-likelihood analysis 
using as fit parameters the released binding energy and the average neutrino en- 
ergy (E Pe ), we find as previous analyses that v e <-> v^ T oscillations with large mixing 
angles have lower best-fit values for (E Pe ) than small-mixing angle (SMA) oscilla- 
tions. Moreover, the inferred value of (Ep e ) is already in the SMA case lower than 
those found in simulations. This apparent conflict has been interpreted as evidence 
against the large mixing oscillation solutions to the solar neutrino problem. In or- 
der to quantify the degree to which the experimental data favour the SMA over the 
large mixing solutions we use their likelihood ratios as well as a Kolmogorov-Smirnov 
test. We find within the range of SN parameters predicted by simulations regions 
in which the LMA-MSW solution is either only marginally disfavoured or favoured 
compared to the SMA-MSW solution. We conclude therefore that the LMA-MSW 
solution is not in conflict with the current understanding of SN physics. In contrast, 
the vacuum oscillation and the LOW solutions to the solar neutrino problem can be 
excluded at the 4<r level for most of the SN parameter ranges found in simulations. 
Only a marginal region with low values of (E Pc ), (Ep^ T ) and E^ is left over, in which 
these oscillation solutions can be reconciled with the neutrino signal of SN 1987A. 
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1. Introduction 

The impact of neutrino properties beyond the standard model on the observed neu- 
trino signal J!] of supernova (SN) SN 1987A has been discussed in numerous works 
during the last 13 years. Most of the early works discussed either mass limits for 
v e or the impact of matter-induced oscillations on the prompt v e burst @. In the 
latter case the MSW effect in the SN envelope would have rendered the prompt v e 
burst unobservable for a large area of mixing parameters assuming a normal mass 
hierarchy. On the other hand, the main neutrino signal of SN 1987A which is the 
v e pulse detected by the reaction v e + p — > n + e + is then influenced only for large 
mixing angles (LMA) by oscillations. It was argued 0, |J that LMA neutrino os- 
cillations would increase the discrepancy between the observed v e spectrum and the 
one predicted by SN simulations @, ||, 0, 0- Moreover, the theoretical prejudice was 
until recently that lepton mixing is, similar to quark mixing, small. Therefore, the 
small mixing angle SMA-MSW solution 1 to the solar neutrino problem seemed to be 
favoured both by the neutrino signal of SN 1987A and by theoretical arguments. 

Meanwhile, the LMA-MSW solution to the solar neutrino problem has become 
the best oscillation solution to the solar neutrino problem [TO] , while Superkamiokande 



atmospheric neutrino data [TT| strongly indicate the need for — > v T oscillations 
with maximal or nearly maximal mixing |1C|| . Such approximate bi-maximal neutrino 
mixing patterns pi can actually arise either in theoretical models based on unifica- 



tion ideas (13| or in "bottom- up" models, such as that of Ref. [14]. It is therefore 



1 This is, for the purposes of the present SN analysis, equivalent to the no-oscillation hypothesis. 
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SMA 


LMA 


VO 


Am 2 /eV 2 
tan 2 $ 


5.0 x 1(T 6 
0.00058 


3.2 x 10~ 5 
0.33 


1.0 x lO" 10 
0.52 



Table 1: Best-fit values for the mass difference squared Am 2 and the mixing angle for 
different oscillation solutions to the solar neutrino problem, from Ref. pQ]. 



interesting to reconsider how serious the discrepancy between the observed and the 
predicted v e spectra of SN 1987A should be taken. 

The main goal of this work is not to obtain those values of the neutrino oscillation 
parameters, tan 2 •§ and Am 2 , that fit best the neutrino signal from SN 1987A. Instead 
we use the best-fit points found from solar neutrino data for three [SMA-MSW, LMA- 
MSW, and vacuum oscillations (VO)] of the four allowed solutions as input and try 
then to assess to what extend the LMA solutions are (dis-) favoured compared to 
SMA assuming certain astrophysical parameters. The case of the LOW solution 
which extends as only one also into the so-called dark side ($ > vr/4) of the solar 
neutrino problem |H| will be discussed separately. There, we are mainly interested in 
the question if the dark side is compatible with the neutrino signal from SN 1987A [[| 
16fl . In all cases, we use as measure for the degree to which the experimental data 
favour one of the allowed solutions over the others their likelihood ratios as well as 
a Kolmogorov-Smirnov test. 

Our main result is that there is a region in the space of astrophysical parameters 
(average neutrino energies (Ei) and released binding energy E^,) predicted by SN 
simulations in which the LMA-MSW solution is either only marginally disfavoured 
or favoured compared to the SMA-MSW solution. For the range of parameters 
compiled in Ref. M, 



14 MeV < (E P .) < 17 MeV 
24 MeV < (E PftiT ) < 27 MeV 
1.5 x 10 53 erg < E h < 4.5 x 10 53 



erg, 



(1.1) 
(1.2) 
(1.3) 



the probability that the LMA-MSW hypothesis is compatible with the detected neu- 
trino signal is extremely varying. We find that for large values of (E Pc ), (Ep ) and 
£?b all three large mixing solutions solutions are practically excluded when compared 
to the SMA-MSW solution. However the probability that the LMA-MSW hypothesis 
is compatible with the neutrino signal rises above 10% for values of {E 9e ), (En) and 
E^ near the lower end of (Ll - L^ ). Moreover, it was recently argued that the average 
neutrino energies are smaller than earlier believed. For example, nucleon-nucleon 
bremsstrahlung that has not been included into the earlier supernova codes tends to 
soften the spectrum ||, a trend which has been confirmed by a recent simula- 
tion 11. This simulation found also for (E 9c ) a rather low value, (E Pc ) « 12 MeV. 
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Since this simulation has not aimed at a self-consistent treatment of the neutrino 
spectra, the values found therein should be taken only as an indication. If (E Pc ) is 
indeed as low as found in Ref. ||, the LMA-MSW solution can be favoured compared 
to SMA-MSW. 

In contrast, the VO and the LOW solutions to the solar neutrino problem can be 
excluded at the 4<r level for most of the range of SN parameters found in simulations, 
and low neutrino energies, {E Pc ) < 12 MeV and (E v ) < 17 MeV, combined with 
low are required to reconcile these oscillation solutions with the neutrino signal 
of SN 1987A. 



2. Neutrino fluences and oscillations 

Massive stars with mass M > 8 M & end their lives in spectacular type-II supernova 
outbursts, releasing almost all their gravitational binding energy via neutrino emis- 
sion [0. Numerical simulations || |6|, |7| as well as analytic considerations Jl8| show 



that this energy is approximately equipartioned between the three neutrino flavours. 
The instantaneous neutrino spectra found in simulations are generally pinched, i.e. 
their low- and high-energy parts are suppressed relative to a Maxwell-Boltzmann 
distribution. The pinching of the instantaneous neutrino spectra is however com- 
pensated by the superposition of different spectra with decreasing temperatures. 
Moreover, we have found that the likelihood function depends only weakly on how 
strongly the spectra are pinched. We present therefore only results using Maxwell- 
Boltzmann distributions for the time-averaged spectra of the neutrinos. Finally, we 
assume that the SN emits the same amount of energy in all neutrino flavours. 

We perform our analysis in the framework of two-neutrino oscillations i> e «-> = 
cos 0^ + sin 0z/ T , which is motivated both by detailed fits of the atmospheric neutrino 
anomaly ||10|| , but also by the results of the Chooz experiment ||19|| . Since the energy 



spectra of z/^ and v T are identical (up to 0(E A G 2 F ) effects), our results do not depend 
on (j). The probability of a v e to arrive at the surface of the Earth can be written as 
an incoherent sum of probabilities, 

P- e - e = P§yPS + PgPi = (1 - P c ) cos 2 d + P c sin 2 d , (2.1) 

where Pj denotes the probability that a v e leaves the star as mass eigenstate Vi and 
Pjf the probability that is detected as v e . The interference terms can be safely 
neglected for all Am 2 of interest because of the spread of the neutrino wave packets 
on their way to the Earth. Deviations from an adiabatic evolution of the neutrino 
states in the SN envelope are characterised by the crossing probability P c . We use 
as approximation for P c the expression valid for an exponential density profile, 

-7Cos 2 i? _ -7 

Pc = — , (2.2) 

l_e-7 K J 
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where the adiabaticity parameter 



7= ^ M1 „ 1=^,1^) (,3) 




is the ratio of density scale height l p and neutrino oscillation length / osc in vacuum. 
The expression (|2.2| ) has in contrast to the one for a linear density profile the correct 
non-adiabatic limit, P c — > sin 2 ^, and describes also for 7 > quite accurately the 
behaviour of P c for a profile typical for a SN envelope, N e oc r~ 3 ||20|| . Furthermore, 



for the case of nearly maximal mixing in which we are especially interested P ss = 
1/2 + 0(tf 2 ). 

The neutrinos from SN 1987A had to cross the mantle of the Earth before they 
reached the Kam and 1MB detectors. Possible matter effects in the Earth can be 
approximated by a box potential, i.e. by replacing PS = |(pj|p e )| 2 with || 

cos 2 -> cos 2 + sin 20' sin(2tf - 20') sin 2 (vrd/Z; sc ) (2.4) 
sin — > sin — sin 2$' sin(20 - 2tf') sin 2 (7rd/C c ) , (2.5) 

where primed quantities are evaluated in matter. The distance d traveled through the 
mantle by the neutrinos and the average density p are different for the Kamiokande 
and 1MB detectors. For Kamiokande we use d = 3900 km and p = 3.4 g/cm 3 and 
for 1MB d = 8400 km and p = 4.6 g/cm 3 @. 

Let us now briefly consider the two extreme limits of the general expressions 
Eqs. ( |2.1| - |2.5|) . For small enough Am 2 , neutrino states evolve in the SN envelope 
strongly non-adiabatically, 7 — > and P c — ► sin 2 Moreover the matter effects in 
the Earth can be neglected and thus the survival probability becomes identical to 
the vacuum oscillation probability, P g g - ► cos 4 $ + sin 4, #. In the opposite limit, the 
evolution of the neutrino states in the SN envelope is completely adiabatic, 7 — > 00 
and P c — ► 0. Thus 

P ss = P^Pfe = cos 2 i? + sin 20' sin(2tf - 20') sin 2 (7rrf/4 c ) . (2.6) 

Hence, Eqs. ( p.l[ - p.2|) describe correctly the survival probability of v e neutrinos in the 
region of parameters favoured by the LMA-MSW and the VO solution to the solar 
neutrino problem and interpolate smoothly for intermediate values of Am 2 and d. 
The v e fluence arriving at the detectors is then 

F Pe = P-eeFl + (1 - P- eE )F» h (2.7) 

where F° stands for the time-integrated flux of neutrinos emitted by the SN. 

3. Likelihood analysis 

The maximum-likelihood method is a particularly well-suited tool for problems like 
the one at hand, where we want to extract the maximal possible information from 
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only 19 neutrino events. We test with the likelihood function 

/ f \ ^obs 

C{a) oc exp I — / n(E, a) dE ) Y\ n {Ei-, a) (3.1) 

the hypothesis that a prescribed neutrino fluence F Pe [E Ve , a) leads to the observed 
experimental data E^ with probability distribution n(E,a). The maximisation of 
C(a) gives an estimate of the values a* = {tan 2 Am 2 , E^, {E Pc ), (Ep b ), . . .} which 
best represent the data set Ei. The confidence region around a* which contains the 
true value of a with a specified probability (3 for k fit parameters is given by 

lii£(a„)-lii£(a) < l - x \ k , (3.2) 

if C{ol) is assumed to be Gaussian near its extrema. 

Apart from parameter estimation given a specific hypothesis, the likelihood anal- 
ysis can also be used to decide which hypothesis fits better the experimental data. In 
this case, the ratio of the likelihood functions for the different hypotheses is a useful 
estimator. We will present the best-fit values for a given oscillation hypothesis in 



Sec. 3.1, while we test which hypothesis fits better the data given a certain range of 



parameters in Sec. 3.2 



The connection between the observed data set Ei and the emitted neutrino flu- 
ence F Pe (Ep e ,a) has been described already in detail in the literature. We follow 
here closely Ref. ||, but use for the detector efficiencies the fit functions given by 
Burrows ||. 

3.1 Best-fit values and confidence regions 

In this subsection we review the best-fit values and confidence regions obtained from 
our likelihood analysis for the different oscillation hypotheses. Since the likelihood 
function we use is identical to the one of Ref. f4|, our results differ only slightly 
owing to small shifts of the best-fit values of the solar neutrino data indicated by the 
most recent analysis |I0|| . Therefore we present only briefly the main points of this 



analysis. 

In the following, the region of SN parameters ( |l.l| - |1.3j ) given in Ref. J7| corre- 
sponds always to the cross hatched region in the figures, while the lower values of 
(E 9e ) found in Ref. are shown as hatched region. The best-fit values obtained for 
the three different oscillation hypotheses are summarized in Table 

The presence of v e «-> z/ MjT neutrino oscillations can affect the average energy of 
v e . We have parameterised the energy as (E Ph ) = t(E Pc ). In Fig. [I], we show 
the 95.4% CL. likelihood contour for the LMA-MSW oscillation hypothesis in the 
E\ 3 -(E i i e ) plane. The neutrino oscillation parameters are given in Table [I] and the 
average energy of is specified by the parameter r. The case r = 1 corresponds 
to the Standard Model or SMA-MSW oscillations, because in both cases oscillations 
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can be neglected altogether. The best-fit point is shifted to smaller values of (E 9b ) for 
increasing values of r, and the 95.4% CL. likelihood contour includes a large portion 
of the hatched region for r = 1.4 and 1, but only touches it for r = 1.7. Parts 
of the cross-hatched region (l.l-i~3|) are however included by none of the oscillation 
hypotheses, and only in the case of the SMA-MSW solution the 95.4% CL. likelihood 
contour touches nearly the boundary of ( |1.1| - |1.3| ). 




Figure 1: Contours of 95.4% CL. likelihood for the LMA-MSW oscillation hypothesis 
as function of E\, and {Ep c ) for r = 1, 1.4, 1.7 and 2. The best-fit points are marked by 
crosses. 

In Fig. |2|, we show the same plot for the VO hypothesis. The oscillation pa- 
rameters are again given in Table |l|. The shift of the best-fit point towards smaller 
values of (E Pc ) is now more pronounced and already for r = 1.4 the 95.4% CL. 
likelihood contour does not include (E Pc ) « 12 MeV. If the VO solution is really 
realized in Nature, then the SN 1987A data would indeed strongly favour smaller 
neutrino temperatures as it is usually assumed. 

The compatibility of the two data sets of Kamiokande and 1MB has been already 
discussed extensively in the literature. Recently, it has been speculated that the 
LMA-MSW solution can improve the agreement between the two experiments |HH . 
We have therefore analysed the two data sets for the SMA-MSW and the LMA- 
MSW case also separately. The separate contours of constant likelihood of the Kam 
and 1MB data sets are shown in Fig. ^| for SMA and LMA-MSW oscillations with 
t = 1.4 and r = 1.7, respectively. As an intuitive measure for the agreement of 
the two experiments we give in Table [| the absolute size of the intersection of fixed 
likelihood regions. If one allows for arbitrary values of E^, the overlapping area of 
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Figure 2: Contours of 95.4% C.L. likelihood for the VO oscillation hypothesis as function 
of E\y and (Ep e ) for r = 1, 1.4, 1.7 and 2. The best-fit points are marked by crosses. 





Kam 


1MB 




oint 




(Eu.) 


E h 


(E* e ) 






E h 


Ax 2 


SMA 


7.8 


5.3 


11.3 


4.8 


10.9 


3.4 


1.9 


LMA, t = 1.4 


6.6 


6.6 


9.9 


5.8 


9.8 


3.6 


2.0 


LMA, r = 1.7 


5.6 


8.8 


7.2 


16.2 


8.5 


4.5 


1.6 


LMA, r = 2 


4.8 


10.8 


5.9 


23.6 


7.1 


5.9 


0.8 


VO, r = 1.4 


6.0 


6.6 


8.4 


7.3 


8.6 


3.9 


1.2 


VO, r = 1.7 


5.1 


7.3 


6.9 


8.3 


7.1 


4.5 


0.4 


VO, r = 2 


4.2 


8.7 


5.7 


9.9 


5.9 


5.3 






Table 2: Best fit values for (Ep e ) /MeV and Sb/lO 53 erg; both for a separate and a joint 
analysis of the two experiments. The difference Ax 2 of the local minima to the global 
minimum (VO solution with r = 2) is also given. 

both the 90% and the 95% C.L. regions increases considerably for increasing r, as 
can be seen from Table |[ Since however the likelihood contours of the individual 
experiments extend up to unrealistic large values of E^, we present results not only 
for the case that E^ can float arbitrarily but also for the case that E^ is restricted to 
lie in the band 1.5 x 10 53 erg < E^ < 4.5 x 10 53 erg. Then, the absolute size of the 
intersection of the two experiments diminishes slightly for increasing r. The same is 
true if one uses the relative size denned as fraction of the intersection with the union 
of the areas. We conclude therefore that the LMA-MSW solution with the current 
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£b restricted 


E h free 




90% C.L. 


95.4% C.L. 


90% C.L. 


95.4% C.L. 


SMA 


6.5 


15 


8.5 


22 


LMA, r = 1.4 


5.7 


14 


8.8 


24 


LMA, r = 1.7 


5.1 


13 


12 


31 


LMA, r = 2 


4.4 


12 


19 


45 



Table 3: Absolute size (in arbitrary units) of the overlapping part of the 90% and 95.4% 
confidence regions of the two data sets. 



best-fit values for Am 2 and $ does not improve significantly the compatibility of the 
two data sets. The compatibility might improve however for a judicious tuning of 
Am 2 and d [|2T] . Finally, we want to comment on the very large values of E^ found 



to be compatible with the 1MB data set (cf. Fig. |3|). Both the energy dependence of 
the cross-section a{v e p — > ne + ) oc i? 2 and of the detector efficiency made the 1MB 
detector blind for v e with E v < 20 MeV and then increasingly sensitive to v e with 
higher energies. Therefore, a decrease in (Ep c ) can be compensated by an increase 
in i?b such that enough neutrinos from the Boltzmann suppressed high-energy tail 
of the distribution are still detected, while the main part of the emitted neutrino 
signal is not seen. Thus and (E Pc ) are, especially in the 1MB data set, strongly 
correlated and the shape of the likelihood contours is very distorted. 

In the analysis outlined above we simply overlaid the theoretical expectations 
for Eb and (E Pc ) with the confidence regions obtained from the experimental data 
set. In the derivation of these confidence regions we have not used any a priori 
knowledge of E^ and (E Pc ). Within the Bayesian approach [p2|, we could combine 



the theoretical prior information p pr ior from Eqs. ( |1.1[ - |1.3|) with the experimental 
information contained in the likelihood function £, 

Pposterior(a) °C £(")Pprior(tt) ■ ( 3 -3) 

If we are, for instance, absolutely sure that the theoretical range B defined in 
Eqs. ( |l.lj - |1.3| ) is correct and assume a flat probability distribution p pr i r inside this 
range for E h and (E Pc ), we obtain 

. . f CCia) for a e B 

<"> = { for (3 ' 4) 

This simple example shows that if one takes serious some theoretical predictions B, 
then it is only important how well the different oscillations scenarios fit the data 
inside B. On the other hand, the behaviour of the likelihood functions outside B 
does not influence at all the posterior probability distribution for a. Furthermore, 
we have noted above that the two fit variables used, (E Pc ) and E^, are strongly 
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A X % 2 = -2 \n(R) 





2 


4 


6 


10 


PhMA = 1 — 0SMA 


50% 


27% 


12% 


5% 


0.7% 



Table 4: Probability j3 that one of the two hypotheses describes the data. 

correlated: Using for instance only the 1MB data set it is impossible to extract any 
useful information about E h . This suggests also to use the a priori knowledge about 
the possible range of E h and (E Pe ) as input in our analysis. 

Another point that can lead to misinterpretations is that the confidence regions 
in Figs. |l] and |2| were constructed relative to the local minima. Therefore, from 
Eq. (|3.2j ), it follows that the large mixing oscillation solutions are relatively penalized 
with respect to the SMA-MSW solution because they have lower x 2 minima. 

3.2 Testing the different oscillation hypotheses 

The ratio of the likelihood functions for different hypotheses measures the degree to 
which the experimental data favour one hypothesis over the other. In order to decide 
how strong the LMA solutions are (dis-) favoured against the SMA-MSW solution 
for a certain range of astrophysical parameters (E^, (E Pc ),t) one has therefore to 
consider the ratio 

R{E h , {E Pc ),t) = — - , (3.5) 

where we treat r as a fixed parameter (r = 1.4, 1.7 or 2 as above). Then the 
probability /3lma that the LMA and not the SMA hypothesis is compatible with the 
data for certain astrophysical parameters is given for f3 — > 1 or (3 — > approximately 
by the usual Ax|. 2 = —2 ln(i?). If both hypotheses have roughly the same probability, 

i.e. Ax% 2 ~ °) tlie equation /?lma(Xlma) + /%ma(Xlma + c ) = 1 nas to be solved for 
a given contour A% 2 = c, cf. Table £|. 

In Eq. ( |3.5D , we compare the likelihood of different oscillation hypotheses for the 
same astrophysical parameters. Therefore, we can now answer questions like "how 
large is the maximal probability that the LMA-MSW solution is compatible with the 
neutrino signal, assuming that the predictions for E h , (E Pe ) and (E„ h ) of a certain 
simulation are correct?" 

In Figs. f| and we show the likelihood ratios \n(R) of the SMA-MSW hypoth- 
esis compared to LMA-MSW and the VO, respectively. Let us discuss first the 
LMA-MSW solution. Restricting the possible E h , (E Pc ), t values to the theoretically 
favoured region ( |1.1[ - |1.3| ), the SMA-MSW solution is always favoured. In particular, 
the LMA-MSW hypothesis can be clearly excluded for large values of E^ together 
with large r. The preference for the SMA-MSW solution is however weak for low val- 
ues of E h , {E Dc ) and r and looses any statistical significance for £ b = 1.5 x 10 53 erg, 
(E 9e ) = 14 MeV and r = 1.4. For (E Pc ) w 12 MeV, e.g. neutrino energies predicted 
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by recent simulations, there are even points in the (EpJ-E^ plane for which LMA- 
MSW is favoured compared to the SMA-MSW solution. Note also that this region 
of relatively small (E Pc ) is favoured by the neutrino signal of SN 1987A for all three 
oscillation solutions. 

The case is different for the VO solution. Within the region ( |1.1[ - |1.3| ), the like- 



lihood ratio is —150 < 21n(i?) < —10. In particular, large values of r and are 
clearly excluded by the data. Only if one allows for rather small neutrino energies, 
(E Pe ) < 12 MeV and r = 1.4, some parameter space remains for which the VO 
solutions is not clearly disfavoured by the data. 

Finally we want to discuss the case of the LOW solution to the solar neutrino 



problem. As has been shown in Ref. [T(| this solution also provides a good global 



fit of the solar neutrino data and extends continuously from the "normal" d < 7r/4 
(light side) case into the region with d > it /4, the so-called dark-side. In order to 
determine the impact of v e <-> u^ tT neutrino oscillations in the LOW region on the 
observed v e signal of supernova SN 1987A we have repeated a likelihood analysis for 
d and Am 2 for this case, choosing for E^, (E Pc ) and r some representative values. In 
Figs. [|and [7|we show the likelihood ratio ln(R) as function of tan 2 $ and Am 2 relative 
to the SMA hypothesis for two values of the average v e energy, {E Pe ) = 12 MeV and 
(Ep e ) = 14 MeV, respectively. The contours of constant likelihood shown correspond 
to ln(i?) = —1, —2, —3, —5, —10, —15, —20, —30, if not otherwise indicated. We find 
that the contours are symmetric with respect to $ = tt/A below 10~ 7 -10~ 8 eV 2 . Thus 
matter effects do not play a role in neither the SN envelope nor the Earth for most 
of the LOW region. Therefore, the dark side is as compatible with the SN 1987A 
neutrino signal as the light side. The absence of Earth matter effects below 10~ 5 eV 2 
also shows that regeneration effects are negligible for the LOW and VO solutions, 
and "explains" why their status is worse than that of the LMA solution from the 
point of view of the present supernova analysis. 

In Figs. |6] and |7| we present also the 90, 95 and 99% CL. contours of the other 



possible solutions to the solar neutrino problem (from Ref. ]T0[). These figures sum- 
marize therefore our findings above. The Earth matter effect regenerate partially the 
v e fluence for Am 2 ~ 10 -5 eV 2 and thereby increases the likelihood of the LMA-MSW 
solution when compared to the LOW and VO solutions. The LMA-MSW solution 
has a large probability to be compatible with the SN1987 data for a considerable 
part of the space of astrophysical parameters, while the LOW and VO solutions are 
generally strongly disfavoured by the data, except for the extreme case low E^, (E 9e ) 
and t. 

4. Kolmogorov-Smirnov test 

We use the one-dimensional Kolmogorov-Smirnov (KS) test as described e.g. in 
Ref. P3[ to estimate the probability that the observed data points E^ agree with 
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the spectral shape predicted by a given distribution n Pe (E). The one-dimensional 
form of the KS test we apply has the disadvantage of being not sensitive to viz. 
the amplitude of the neutrino signal. This disadvantage is however compensated 
by its property to be distribution-free that makes the one- dimensional KS test at- 
tractive. In Fig. H we show the probability that the signal observed in the two 
detectors is consistent with the SMA-MSW hypothesis, as function of (E Pe ). Per- 
forming this test for the two experiments separately, we find a small region of overlap 
at (E 9e ) ~ 11 MeV where this hypothesis has a probability of more than 10% in both 
experiments. The combined data set has a probability higher than 10% in the broad 
range (E Pc ) = 9 MeV - 14 MeV, with a maximal probability of 46%. Since the 
maximum occurs at the same energies as the best-fit values found in the maximum- 
likelihood analysis, we conclude that the essential information is contained in the 
form of the spectra, not in the absolute number of events. 

In Fig. H we show the same test for the LMA-MSW hypothesis. Its probability 
as function of {E Pc ) has the same shape as in the case of SMA-MSW but is shifted 
to lower energies by 0.7 MeV (r = 1.4), 1.7 MeV (r = 1.7) and 2.8 MeV (r = 2), 
respectively. Consequently, an exclusion of the LMA-MSW solution because of its 
preference of low (E 9e ) requires a correspondingly precise knowledge of the average 
neutrino energies emitted by the SN. The maximum of the probability as function of 
(E Pe ) increases from 50% (r = 1.4) to 65% (r = 1.7) and 80% (r = 2). Moreover, we 
find that the LMA-MSW hypothesis is consistent with the data at the 5% (r = 1.4), 
2% level (r = 1.7) and 1% level (r = 2) for (E Pa ) = 14 MeV, while the compatibility 
increases to 25% (r = 1.4), 16% level (r = 1.7) and 5% level (r = 2) at the lower 
end of predicted neutrino energies, (E Pc ) = 12 MeV. 

Finally, we present in Fig. [10] the results of the KS test for VO. Now the maxima 
of the probability as function of (E Pc ) compared to the SMA-MSW solution are 
shifted by 2.2 MeV (r = 1.4), 3.7 MeV (r = 1.7) and 5 MeV (r = 2), respectively. 
The maximum of the probability increases from 59% (r = 1.4) to 68% (r = 1.7) 
and 70% (r = 2). Thus the most probable value of {E Pe ) for r = 2 is reduced by a 
factor two compared to SMA-MSW and also the cases r = 1.7 and r = 1.4 are highly 
incompatible with the range of neutrino temperatures predicted by simulations. On 
the other hand, VO together with rather low neutrino temperatures, (E Pc ) « 12 MeV 
and r = 1.4 are at the 2a level compatible with the spectral shape of SN 1987A signal. 



5. Discussion and Conclusions 

Recently, several papers discuss the influence of neutrino oscillations on the signal of 
SN 1987A. Reference |H| stresses that matter effects in the Earth partially regenerate 



the v e flux in the case of LMA-MSW oscillations. This effect can be clearly seen in the 
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plots of R in the Am 2 -tan 2 ■& plane 2 , Figs. |-|7], where the likelihood of the LMA-MSW 
solution increases for Am 2 ~ 10~ 5 eV 2 . Furthermore, Ref. |ZT| finds that the Kamll 
and 1MB data sets become more compatible to each other for LMA-MSW oscillations 
because of the different Earth matter effects seen in the two detectors. The general 
tendency of their findings agrees with ours although there are some differences in 
the details. A possible explanation for the differences is the different approach used: 
Ref. [21 1 binned the 1MB data into two bins while we used a continuous likelihood- 
function together with an experimental energy resolution function. Reference [f24 | 
concentrates on the question if an inverted mass hierarchy m„ 3 < rn ui < m U2 with 
Am 2 2 = AtjIq and Am 2 3 = Am 2 tm can be excluded by the SN1987 signal. They 
conclude that the SN 1987A signal disfavours this scheme unless \Ui^\ 2 < a few 10~ 4 . 
Finally, the authors of Ref. p5 consider the compatibility of the LSND hint for 
<-> v e oscillations with the SN 1987A signal and find that the LSND result is 
disfavoured. 

In this work, we have put the main emphasis on scrutinizing whether or not 
and to what extent large mixing oscillation solutions to the solar neutrino problem 
are really excluded by the SN 1987A signal. Our main result is that the LMA- 
MSW solution is not in conflict with the current understanding of supernova physics. 
In contrast, the VO and the LOW solutions to the solar neutrino problem can be 
excluded at the 4<r level, for most of the range of SN parameters found in simulations. 
Only a marginal region with low values of (E Pc ), (E v ) and is left over, in which 
these oscillation solutions can be reconciled with the neutrino signal of SN 1987A. If 
either of the large mixing solutions will be finally established by future solar neutrino 
experiments, the neutrino signal of SN 1987A will allow one to restrict severely 
the possible range of the neutrino energies and the released binding energy in SN 
explosions. 
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Figure 3: Contours of constant likelihood for the SMA-MSW (top) and the LMA- 
MSW oscillation hypothesis as function of and (E Pc ) for r = 1.4 (middle) and r = 1.7 
(bottom). The two experiments are analysed separately. 
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Figure 4: Likelihood ratio ln(i?) of the LMA and SMA oscillation hypothesis as function 
of E h and (E Pc ) for r = 1.4 (top), r = 1.7 (middle) and r = 2 (bottom). 
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Figure 5: Likelihood ratio ln(R) of the VO and SMA oscillation hypothesis as function 
of E h and (E Pc ) for r = 1.4 (top), r = 1.7 (middle) and r = 2 (bottom). 
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Figure 6: Likelihood ratio ln(i?) relative to the SMA-MSW hypothesis, as function of 
tan 2 i? and Am 2 /eV 2 for r = 1.4 (top), r = 1.7 (middle) and r = 2 (bottom) together with 
the 90, 95 and 99% CL contours for the different solutions to the solar neutrino problem 
from Ref. 0. All figures for E h = 1.5 x 10 53 erg and (E Pc ) = 12 MeV. 
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Figure 7: Likelihood ratio ln(R) relative to the SMA-MSW hypothesis, as function of 
tan 2 i? and Am 2 /eV 2 for r = 1.4 (top), r = 1.7 (middle) and r = 2 (bottom) together with 
the 90, 95 and 99% CL contours for the different solutions to the solar neutrino problem 
from Ref. 0. All figures for E h = 3 x 10 53 erg and (E Pc ) = 14 MeV. 
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Figure 8: Probability p as function of (Ep c ) given by the KS test that the observed data 
by Kam, 1MB and the combined data set are consistent with the SMA-MSW oscillation 
hypothesis. 
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Figure 9: Probability p as function of (E Pc ) given by the KS test that the combined data 
set is consistent with the LMA-MSW oscillation hypothesis for r = 1.4, r = 1.7 and r = 2. 
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Figure 10: Probability p as function of {Ep c ) given by the KS test that the combined 
data set is consistent with the VO hypothesis for r = 1.4, r = 1.7 and r = 2. 
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